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Abstract 

The confined, quasi-two-dimensional guiding center plasma and a system of interacting line vor- 
tices in an ideal fluid are examples of Hamiltonian systems with infinite interaction distances. The 
existence of metastable states with negative specific is investigated by standard entropy maximiza- 
tion of the thermodynamic limit of vortices as they become infinitesimal and form a continuous 
field. We find metastable states and suggest that these imply a runaway reaction leading to a 
rapid expansion of a confined plasma or fluid similar to the rapid collapse of globular clusters in 
astrophysics. 
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I. INTRODUCTION 

There has been considerable interest recently in metastable states and negative specific 
heat, particularly related to fluids and plasmas [4, 20]. The guiding center plasma or ideal 
fluid vorticity model for quasi-2D columns of electrons or lines of vorticity is similar to 
the widely studied two-dimensional model but the lines contain small variations which can 
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change the dynamics of the system. An ensemble of filaments {0i(t), . . . , 0jv(t)}, 



where 1/e is the coupling constant, all the filaments have the same average circulation, T, a 
is a core elasticity constant related to the frequency, and / is the length of the period under 



in the plane with a parameter r representing the third dimension. This is under special 
asymptotic assumptions that the filaments are nearly parallel and far enough apart [16]. 
Because of the nearly parallel assumption, this model neglects vortex stretching which is 
assumed to be too small to affect the statistics. An example of such an ensemble is pictured 
in Figure 1. 

There is a significant difference between the two-dimensional one-component Coulomb 
plasma and this quasi-2D model. Because the total energy is entirely dependent on how far 
apart the lines are, an ensemble of two-dimensional lines at a fixed energy has a maximum 
radius beyond which the lines cannot move while the quasi-2D lines can move, theoretically, 
as far apart as there is space available because the potential energy can always be balanced 
by altering the kinetic energy. In addition, we show in this paper that metastable states 
exist in the quasi-2D model that do not exist in the 2D model. 

Statistical derivations for three dimensional fluids of any kind have tended to focus on sin- 
gle vorticity columns [6, 7] and [8, 9] or statistical treatments of ensembles of two-dimensional 
point vortices [1, 10-12]. The nearly parallel vortex filament model of [15, 16] for Navier- 
Stokes fluids is an exception. 

The nearly parallel model can be extended to electron columns by a generalized vorticity 
model for electron plasmas [5, 13, 14] which takes the magnetic and electric fields into account 
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periodic boundary conditions 0«(O) = (f>i(l). Each filament 4>i(j) = (rc^r), yj(Y)) is a vector 
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FIG. 1: This output from a Monte Carlo simulation of the energy functional. 

as well as the vorticity. In the case of charged particles, the vorticity must, essentially, be 
gauge invariant: the magnetic field, —eB/m = (— e/m) V x A, and the charged fluid vorticity, 
u = V x v , combine into a general vorticity field Q = m _1 V x p where the generalised or 
"canonical" gauge invariant momentum is, p = mv — eA, m is the electron mass, — e is the 
electron charge, v is the fluid velocity field, and A is the magnetic vector potential field. 
The generalized angular momentum is L = r x p. 

Electron column core sizes are small, equal to the Larmor radius, a, about a tenth of a 
millimeter for electrons, hence the core structure may be abstracted by the local induction 
approximation (LIA). The size of the local region, which is the wavelength of the highest 
energy frequencies on the filaments, is naturally given by the London wavelength, h = c/u pe , 
where the electron plasma frequency is given by ut e = ne 2 /meo and n is the electron number 
density of the individual filament [13] . Because the filaments are all nearly parallel, collisions 
between points not in the same cross-section are small, and the two dimensional Coulomb 
interaction gives the potential to leading order. 

II. METHOD 

A filament nearly parallel to the z axis has a C 2 curve 0(r) = (x(r),y(r)) with L 2 ^ 
derivative where / 1 is a small length scale. Assuming no short-wave disturbances (knots, 
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kinks, etc.), the filament has a LIA kinetic energy functional, 

where T is the generalized circulation of the filament [13] and a = \og(b/a) + 1, b is the 
arclength of the "local" region where the induction takes place, i.e. the wavelength of the 
highest frequency modes on the filament, and a is the core size of the filament. The ratio 
b/a does not vary significantly over the length of the filament; therefore, it can be assumed 
to be constant [9]. The core size a of electron filaments is the Larmor radius while a natural 
choice for b is the London wavelength, b = c/u pe , where the electron plasma frequency is 
given by u 2 e = ne 2 /me and n is the electron number density of the individual filament [13]. 

For N such filaments, they have, under suitable asymptotic assumptions such that the 
filament core size is much smaller than the intervortex spacing, a 2-D Coulomb interaction 
such that, 

2 

-I^I*log|&-&-|, (3) 

i<j 

where 1/e is the coupling constant and all the filaments have the same average circulation 
[14]. The kinetic energy is the square amplitude of the filament with zero amplitude im- 
plying zero kinetic energy; hence the filament's self energy is not directly connected to the 
microscopic temperature. 

The kinetic energy of generalized angular momentum for a single filament of electrons is 
A 1 = \loo 2 where / is the moment of inertia of the filament, / = j l Q dr |0| 2 , and u is the 
generalized frequency of rotation. For N filaments, all with the same frequency, 

A N [$i---M^y j Tdrrl^l 2 , (4) 
1 i J ° 

where fi'T = u 2 . Now choose mass and charge units such that e = m = 1 for electrons. 

No boundary conditions are imposed perpendicular to z since the plasma is fully contained 
by the magnetic field and never contacts any surfaces. The confinement radius R is a 
bounding radius representing the distance from the center of the plasma to the outer edge 
where the density falls to zero (sometimes abruptly). This is smaller than the radius of 
the container. Provided the magnetic surfaces to which they are confined are closed, the 
filaments cannot interact with material surfaces. 




E N [f 1 ...f N ] = Y t f ^ 
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In the thermodynamic limit as iV — > oo such that total circulation, A = TN, is constant, 
the filaments have area density /(c, r, r) = Tg(c, r, r) (where c = dtf/dr and T is the 
generalized circulation) such that, if a = (c, r), the kinetic energy is, 

T=\j dVdr/ac 2 , (5) 

where c = ||c||. In a rotating frame the kinetic energy of angular momentum becomes a 
potential added to the interaction potential, 

V = dVdr/|rf - iy dVdVdr/Z'loglr-r'l, (6) 

where /' = f(c',r',r) and particle number is N = J dVdrg. All integrals are over the 
interior of the torus. Since the energy functional does not depend on r, the density with 
maximal entropy does not depend on r either (as one can show from the variation); therefore, 
we drop the integrals over r and assume the area density is constant in r. All functionals 
are now per unit length. 

For a fixed energy system, the entropy (with Boltzmann's constant ks = 1), 

S = -Jd 4 ag\ogg, (7) 

is maximal in the most-probable macrostate. 

To maximize the entropy, we must solve the variational problem, 

5S = 0, (8) 

subject to the constraints, T + V = E, N fixed. That is for some small parameter A, we 
define a family of density functions g(cr; A), such that entropy is maximal (SS = dS/dX = 0) 
at g(a;0). 

The method of Lagrange multipliers provides the equation for the variation of S subject 
to the constraints, 

5S + p'SE + v5N = 0, (9) 

where /3' is inverse temperature, v is a normalization parameter. (Angular momentum is 
automatically conserved by the rotational invariance of the energy.) Taking the variation 
([3]) we have, 

log g + 1 + pT (|c 2 + |r 2 + ^ + v = 0, (10) 



where 



^(f) = ~J dfp(f) log |f - f\ (11) 



is the 2D Coulomb potential. Let H = f c 2 + ^r 2 + tp. 
Solving 10 for the density, 

g = Ae~^ H , (12) 

where A = exp[— [y + 1)] is a normalization constant that gives the particle number N. Let 
j3 = j3T. Showing the equipartition of energy, the average kinetic energy per filament is 

Jd 2 c/fc 2 _ /d 2 cexp(->c 2 ) |ac 2 _ 1 
Jd 3 c/ Jd 2 cexp(-±/3ac 2 ) 1 ' 

To solve for the spatial density, p, we can integrate equation 12 over all "velocities" , c, 

P = J d 2 c/ = J Bexp[-/3(^ + /x , r 2 /2)], (14) 

where B = A{2tt/(3). 

Replacing p in the potential with the above equation gives an integral equation for the 
most-probable potential inside the circle, 

ij){f) = - J dr'Be-^ )+ »' r2/2) log |f - f\. (15) 

This integral equation is equivalent to the Poisson equation, 

W(0 = ~-P = l -^^-^ + ^ 2)l/£ Ir1 < R (16) 
e [ \f]>R 

with boundary conditions such that ip and dip /dr are continuous at the boundary r = R. Be- 
cause of the axisymmetry of the energy the potential must also be statistically axisymmetric. 
Converting to polar coordinates and 16 gives the ODE, 



1 d 

r dr 



(M\ = -i^e-^WHV^ (17) 
\ dr J e 



for r < R. 

Because the potential is repulsive, the only solutions to this equation have finite den- 
sity everywhere (unlike gravitational systems, infinite densities such as black holes are not 
possible in a repulsive Coulomb system). To simplify, we make the change of variables, 
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vi = /3(ip — 1/>(0)) and r\ = ^4^/3 A exp[— f3ip(0)]/e ■ r = a/ ' Atx j3 p(0) / e ■ r , which simplifies the 
ODE (17) to, 



+ I** + ^-MW/W = (lg) 

where z = ^ An p(ti) (3 / eR and with boundary conditions, 

Vl (0) = v[(0) = 0, (19) 

where = dv\jdr\. Note that the plasma density can be written p = p(0)e~ Vl or {v\ = 
— log[p/p(0)]); therefore, the variable V\ describes how the density changes as the distance 
from the origin changes and decreases monotonically from at r± — 0. 
The viral theorem of Clausius may be applied, 

2T+V = 3pV, (20) 

where 

P= f d 2 c/i«c 2 , (21) 

Jr=R 6 

is the surface "pressure" and V = ttR 2 is the area of the circle. From the equipartition 
theorem, where A is the total circulation, 

r=j (22) 

Using the Virial theorem, 

E = 3pV — T = 3pV — (23) 

and 

V = 3pV - 2j. (24) 

Evaluating the surface pressure, 



where 

/ = Aexp 

Integrating 17, circulation A is given by, 



(5 ( a— + p'—+ij 



(26) 



7 ' 
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Then 

= ezv[(z) 
1 A ' 

For the rest of this section, v[ and v\ shall refer to v[(z) and vi(z) only. 
The pressure is, 



(28) 



_2p{R) _2 p(0)e-^ 2 / 2 
Since r 2 /r 2 = \4n(3p(0)/e\, and eliminating instances of f3 with 28, 

2 e ^2 e -yi-/Vil 2 /2 ^2 2 2 e -i;i+ M '^e_R 2 /(2A) 

P= 34^R2 ^ = Q^We (-^i) 2 ' ^ 



Let fi = e/i'i? 2 /(2A), then 



and the energy from 23 is 



A 2 2 -v 1 +ixzv[ 



A 2 / z ^ e -^+^i i 



The entropy, obtained from 7, is given by, 



5= M/3 — logi?) - A logp/3 2 + A log 2ixT \ . (33) 



r r V e 

The specific heat is, 



,/~,2 / d2 



c = _ = ^ (34) 
dT dT ' 



When fxzv[r / R = V\{z) we have a constant density solution, where the potential energy 
and the confinement are perfectly balanced. For \izv' x r 2 jR 2 > v\(z) the potential dominates, 
and the density tends to favor expansion with higher density toward the wall of the container. 
These profiles, although stable, are not suitable for containment. For fizv[r 2 /R 2 < v\{z) 
the magnetic confinement dominates and a more Gaussian density profile is preferred with 
higher density toward the middle. These profiles are useful for containment. 

Numerically evaluating v\ and v[ for a range of z values (Fig. Ill) in the strong rotation 
regime fi > 0.5, we have negative specific heat. At // = 0.5 the energy is zero indicating 
that the potential energy, which is the only energy component that can take negative values, 
perfectly balances the kinetic energy and the confinement potential; at this point the specific 
heat is zero. In the weak rotation regime with \i < 0.5, the specific heat remains positive 
(Fig. 3). 
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III. RESULTS 



Entropy maximization of the model yields the following results: 

1. An expression for the energy in terms of z, vi(z), and v[(z), 

A 2 ( z i e -w< 1 \ 

E -T\-w^~r^))^ (35) 

where A is the total circulation, e is the coupling constant for vorticity, and ji is a 
parameter determining strength of the angular kinetic energy; 

2. An expression of the inverse temperature of the filaments (Boltzmann's constant, ks = 

i), 

1 =/ , = _£2LM; (36) 

3. An expression for the central density, 

"<°> = 4$?> < 37 > 

where R is the radius of the confinement area; 

4. A second ODE governing the density at a particular energy/temperature, 

^1 + I*!! + e -v l{ri)+ ^ v[{z y z = (3g) 
ar{ T\ dr\ 

with boundary conditions, 

Vl (0)=ui(0) = 0, (39) 

where = dv\jdr\\ 

5. And an expression for the density: 

p(r) = p(0) exp(-ui(r) + fir 2 zv[{z) / R 2 ) , (40) 



where r = r±/ ^ 'Air j3 p(0) / 'e is the distance from the center. 

From numerical evaluation of 38 the energy (35), inverse temperature (36), and core density 
(37) are computed as functions of z. These are plotted for several values of \x (Fig. Ill, Fig. 
3, Fig. 4). 
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FIG. 2: For regimes with /i < 0.5, confinement is weak enough that stable states exist, indicated by 
energy decreasing with z, but beyond that all states are metastable, indicating that equilibrium is 
not possible. When energy increases with z, the system is metastable, and, because this leads the 
system to run away to hotter temperatures where the central density decreases, runaway expansion 
results. 

IV. DISCUSSION 

Negative specific heat indicates a metastable state where the entropy, rather than being 
globally maximal, is at a saddle point or local maximum, and the system evolves out of 
this state by increasing temperature, where the term "temperature" refers to the amount of 
kinetic energy per unit circulation and not the electron temperature. Because the specific 
heat is negative, increasing temperature lowers the total energy. In gravitational systems 
where this phenomenon was first described, this results in gravothermal catastrophe where 
a globular cluster, for example, experiences a core collapse. The collapse of the core results 
in reduced gravitational potential and increased stellar velocities. Thus, the overall energy 
decreases while the temperature increases. It has been argued by analogy that negative 
specific heat in magnetically confined, neutral electron-positron plasmas would result in 
core collapse leading to the possibility of nuclear fusion as columns of electrons and ions 
collapse into one another [4]. In a previous paper, however, we showed that self-energy 
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Temperature vs. z 
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FIG. 3: Temperature decreases with z. Therefore, if energy increases with z, the specific heat is 
negative. If energy decreases with z the system has stable equilibria. 

causes an anomalous expansion in the mean radius of the density profile and developed 
a mean-field formula for it, confirming it with Monte Carlo simulations [19]. The same 
approach in a microcanonical ensemble shows that the system exhibits negative specific 
heat [20] . Therefore, while magnetically confined quasi-2D plasmas exhibit negative specific 
heat indicating metastable energy states where no equilibrium exists, because the potential 
is repulsive, the instability is a runaway expansion. 

This expansion arises because increases in kinetic energy and decreases in potential energy 
both cause expansion. Analyzing Figs. Ill, 3, and 4, the mechanism works as follows: (1) 
The potential energy decreases. If the magnetic confinement is sufficiently weak at the 
given density, expanding it decreases the potential energy. (2) The kinetic energy increases 
in response, but, because the system has expanded, energy is lost at the outer edges of 
the plasma. Thus, the filament temperature has increased but energy decreased. (3) The 
increased kinetic energy causes the density profile to spread (by increasing the variance), 
causing further expansion, and returning the cycle to step (1). As the potential energy 
continues to decrease, the kinetic energy continues to increase, and the system experiences 
a runaway expansion directly analogous to the runaway collapse of a gravitational system. 

Unlike in the gravitational systems, however, because the specific heat is negative for all 
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FIG. 4: As increases, confinement becomes stronger and the central density decreases less with 
increasing z, but, because it is decreasing in the metastable regime, this indicates a core expansion. 

energy values, no equilibrium states exist and the plasma never settles. Because the system 
is forced, like a forced snow-pile (with continuous snow fall), as long as the metastable state 
persists, the energy input into the system from the externally applied currents balances the 
energy loss. When the metastable state ceases, however, the rate of energy loss increases 
dramatically and the rate of energy input no longer balances the rate of energy loss. This 
expansion causes a dramatic loss of core electron temperature which has been observed in 
experiment [21], a direct result of heat transport from the core to the outer edge via the 
runaway expansion mechanism. Recovery occurs when the filaments expand far enough 
apart that the rate and magnitude of their collisions decreases, and the expansion slows. 
Supposing that complete plasma disruption does not occur, the input current allows the 
energy to increase again, returning the plasma to a metastable state. Hence, the instability 
is regular and repeated. 

This instability is a direct result of the confinement as Fig. 3, where confinement is weak, 
shows. In the strong confinement regime, expansion decreases the potential because V\(r) 
increases from the core showing that density (40) decreases; when confinement is weak, 
however, expansion increases it because v\{r) decreases from the core, and any further 
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expansion increases the density at the edge of the confinement area. In a real system, 
weakly confined plasmas do not persist for long before complete disruption however, so this 
regime is not useful for sustained fusion. 

Negative specific heat has been found in neutral, two-component, electron-positron plas- 
mas [4]. This model, however, ignores Coulomb interactions between charged particles and 
does not take into account that the positively charged ion lines move far more slowly than 
the negatively charged electron lines. For a neutral electron/ion plasma, a two-fluid magne- 
tohydrodynamical model is appropriate. At 100 million degrees K (the minimum target of 
tokamak reactors although we are not assuming a toroidal geometry here) electrons travel 
at a mean velocity of 40,000 kps while deuterons travel at "only" 600 kps. Instabilities in 
tokamak plasmas such as the sawtooth instability may be related to the metastable state 
described in this paper [2]. Future research will focus on toroidal geometries to determine 
if this is likely. 
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